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The theory for Unear deformations of fluid microparticles in a laser beam of Gaussian profile is presented, when 
the beam focus is at the particle center as in optical trapping. Three different fluid systems are considered: 
water microdroplet in air, air microbubble in water, and a special oil-emulsion in water system used in exper- 
iments with optical deformation of fiuid interfaces. We compare interface deformations of the three systems 
when illuminated by a wide (compared to particle radius) and narrow laser beams and analyse differences. 
Deformations of droplets are radically different from bubbles under otherwise identical conditions, due to the 
opposite lensing effect (converging and diverging, respectively) of the two; a droplet is deformed far more than 
a bubble, cetera paribus. Optical contrast is found to be of great importance to the shape obtained when 
comparing the relatively low-contrast oil-emulsion system to that of water droplets. We finally analyse the 
dynamics of particle motion when the laser beam is turned on, and compare a static beam to the case of 
a short pulse. The very different surface tension coefficient implies a very different time scale for dynamics: 
microseconds for the water-air interface and tens of milliseconds for the oil-emulsion. Surface oscillations of 
a water microdroplet are found always to be underdamped, while those of the oil-emulsion are overdamped; 
deformations of a microbubble can be either, depending on physical parameters. 
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1. Introduction 

Optofluidics, the marriage of optics and microfiu- 
idics, is a field wliicli has experienced rapid growth 
in recent years. Many areas of apphcation have al- 
ready been demonstrated, including chemo-biological 
applications [ll-H, solar energy applications micro- 
droplet lasing t5|, optically controlled droplet transport 
and coalescence [61 which can be a vehicle for DNA cal- 
culations and a considerable range of dynamically 
configurable integrated fluid-based optophotonic devices 
(8h10I| . Lab-on-a-chip applications of microfluidics have 
already been realised for various purposes d, [llj , For an 
in-context review, see Ref. [l^ . 

Applications of optical trapping and tweezing of mi- 
croscopic fluid droplets are easy to imagine. Optical 
tweezers allow detailed manipulation of microsystems 
[ist . A laser can for example measure the size of a spher- 
ical droplet via its Mie resonances (an old idea, see e.g. 
Ref. [ij) and automatically choose and transport the 
desired quantities of reactants for microchemistry on a 
chip. A related technique was alredy demonstrated for 
sorting of particles by size or refractive index jl5ll] . 

In light of the technological promise of laser manipula- 
tion of microfluidics, it is of obvious interest to consider 
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the deformation and motion of a microdroplet moved by 
a laser. The classic degrees of freedom when manipu- 
lating microobjects with light are pushing and trapping, 
and recently it was discovered that special laser beams 
can also pull against the direction of propagation Tg']. 
While the fluid droplet behaviour under pushing and 
pulling have been analysed in the past [17[, the present 
paper has particular relevance to the case of optically 
trapped droplets, bubbles and emulsions and, even more 
so, deliberate deformation of droplets, bubbles and emul- 
sion by a Gaussian beam. 

As well as using lasers to move fluids around, con- 
trolled deformation of fluids using lasers has attracted 
considerable interest. A number of impressive experi- 
ments in recent times demonstrate the potential for us- 
ing laser light for detailed manipulation of microflows. 
In particular, two-fluid systems of oil emulsion droplets 
have been used, in which the surface tension of the 
liquid-liquid interfaces can be reduced to a millionth of 
that of water-air. References include a series of works 
by Delville's grouplHIll] and associated theory [H,!!!], 
as well as "optical sculpting" work at the Central Laser 
Facility [25l - l27| . Mitani and Sakai used a similar system 
for the measurement of surface tension [28j . 

The manipulation of liquid surfaces with laser light 
goes back to Ashkin an Dziedzic [2^. It was soon con- 
flrmed [sO] that a correct description of the observed de- 
formation could be given by integrating the Minkowski 
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electromagnetic stress tensor across the surface. 

A classical experiment on the deformation of droplets 
by laser pulses was that of Zhang and Chang [3l| , whose 
results could be satisfactorily analysed with essentially 
the same means [32:, , ,33] , drawing on Mie's theory for 
light scattering on spherical particles [13] • 

The deformation of a liquid droplet with ultra-low sur- 
face tension trapped in a laser beam was studied exper- 
imentally by M0ller and Oddershede [HI As the droplet 
is elongated in the direction of the laser propagation, 
its radius in the cross-beam plane decreases. The ra- 
dius was measured as a function of laser power for three 
different droplet sizes and compared to a simple model. 
A linear decrease with laser power, predicted by their 
model as well as the one presented herein, was observed 
for low laser powers, but beyond a threshold power the 
variation with power becomes much weaker. While still 
limited to small deformations, the theory presented in 
this paper is able to elucidate the observations by M0ller 
and Oddershede. In particular, their experiment gives 
an indication of the range of validity of the linear defor- 
mation theory. 

In the present manuscript we also study the "oppo- 
site" system exemplified by air bubbles in water, such 
as considered experimentally in [s^ . The theory for de- 
formation of bubbles is virtually identical to that for 
droplets, but the deformation is qualitatively different. 
In general the case denoted "bubble" includes all cases 
where the illuminated spherical fluid object has a lower 
refractive index than the surrounding medium, as is the 
case for air bubbles in water. In section [2] where we 
lay out the theoretical framework, we shall only speak 
of "droplets" for simplicity, in the understanding that 
the exact same expressions are applicable to the bubble 
situation. 




Fig. 1. The set-up considered. An initially spherical parti- 
cle of fluid 1 enclosed within a fluid 2 is illuminated by a 
Gaussian beam so that the particle's center is at the beam 
focus. 

We consider the set-up shown in Fig. [TJ A droplet 
(bubble) of radius a when unperturbed, density p2 and 
refractive index 712 sits in a medium of density pi and 
refractive index ni and is illuminated by a laser beam 



whose focal waist is wq which can be bigger or smaller 
than the droplet radius. The droplet is deformed by the 
optical force acting upon it. We consider both the case 
where the laser is pulsed, and, with certain approxima- 
tions also the case where it is trapped in a static beam. 

We work in the analytical framework made possible 
by Mie's scattering theory [33]. Analytical results for 
the electric field near a deformed sphere do not exist 
at present, although certain considerations have been 
made, e.g., for the change in lasing properties [stI IbsI]. 
We will work with the internal and scattered electric 
laser field as for a spherical droplet whereas all other 
droplet dynamics take account of deformations to linear 
order. For the case of a short pulse, this is virtually 
exact since hydrodynamic motion takes place after the 
duration of the pulse, while for the case of a trapped 
droplet it is an approximation, the best such presently 
available short of direct numerical simulation for which it 
is harder to keep track of the individual physical effects. 
In the case of oil-emulsion systems the dielectric contrast 
is small and the difference in internal diffraction pattern 
will be modest due to a small deformation, whereas for 
higher contrast (e.g. water in air) a fully extended theory 
is certainly desirable and we plan to return to this in 
future work. 

The theoretical framework of both the fluid mechanics 
and the optics of the problem at hand are laid out in Sec- 
tion O A number of numerical examples are shown and 
discussed in section[3]whereupon concluding remarks are 
made. The effects of gravity (including buoyancy) are 
neglected throughout the manuscript. 

2. Theoretical framework 

We consider the case of an incident field impinging on 
a droplet with radius a and refractive index rti, embed- 
ded in a medium of refractive index rt2. In the present 
paper we shall assume rii and n2 to be real for simplic- 
ity, although the generalisation to a droplet of complex 
rii, hence absorbing droplet, is straightforward. We as- 
sume all fluids to be nonmagnetic, i.e., having vacuum 
permeability. 

2. A. Electromagnetic force density 

In general, the electromagnetic (EM) force density act- 
ing in an isotropic dielectric medium may be written [39| 



1 d 



dt 



{£ X H) 



(1) 



A discussion of the three terms of Eq. ([T]) for the present 
context is found in Ref. ^32] ■ The last term of Eq. ([T]) is 
called the Abraham term. Its interpretation has been 
debated for a century [11], but this is of no concern 
in the present investigation; since the time variations 
of our EM fields are at optical frequencies, immensely 
quicker than the mechanical response time of the fiuids. 
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the Abraham term averages to zero over an optical pe- 
riod and is of no consequence, hence we omit it in the 
foUowing. 

The second term of Eq. ([T]) is called the electrostrictive 
term [3§|. It may be interpreted as a tendency for the 
material medium to be attracted to areas of higher field 
intensity. When the EM field is switched on abruptly, 
the fluid will be slightly compressed in areas of higher 
intensity and vice versa. Soon the material itself will 
counter this compression process with an elastic counter- 
pressure until mechanical equilibrium is reached on a 
timescale of the time a sound wave takes to traverse 
the fluid. Once equilibrium is reached, the electrostric- 
tion term is exactly cancelled by a corresponding rise in 
hydromechanical pressure and is absorbed therein (this 
point was discussed in the present context in Ref. [32]). 
What remains is a very slight non-uniformity of the fluid 
density due to the pressure distribution giving rise, in 
principle, to a slightly nonlinear optical response which 
we shall neglect. Ignoring the strictive term as we shall 
in the following thus means the theory is valid only for 
timescales longer than the traversal time 2a/cs, Cg be- 
ing the speed of sound. For water and air (typefying a 
liquid and a gas) this means that our timescale must is 
restricted to being greater than 



t >1.4ns- 



1/im 



t >5.9ns- 



1/im 



liquid droplet 
gas bubble. 



When considering very short pulses this should be born 
in mind (for further discussion see Ref.^d|). For the 
present endeavour we shall assume all pulses to be much 
longer in duration than these timescales. 

A final note on the EM force density concerns a sug- 
gestion by Peierls [4l| that the gradient term (first term 
of Eq. ID)) should contain a quadratic term in the refrac- 
tive index, oc (e — 1)^. The nonlinear term was expressly 
omitted in [33|, where its effect was estimated to 10% for 
water. What to make of this has been a matter of de- 
bate, however, as it was found to be inconsistent with an 
experiment by Jones and Leslie [i^l, as also detailed in 
the appendix of Ref. [s^ . A force density nonlinear in re- 
fractive index signifies a nonlinear optical effect, hinting 
that the term has to do with the material's compress- 
ibility. Indeed, the apparent mismatch involved would 
appear to have been resolved by Lai et al. [43] who show 
that, just as for electrostriction, the nonlinear force term 
is exactly cancelled by an increase in pressure, as long 
as the latter has time to build. It is a good example of 
how subtle it can be to go from microscopic to macro- 
scopic description in systems where electromagnetic and 
mechanical forces act together. 

We shall therefore consider only the first term of the 
force density H]), often referred to as the gradient force, 
which gives a contribution only at the fluid surface, 
where e varies. It is equivalent (e.g. [ij] §59) to inte- 



grating the gradient of the Maxwell stress tensor 

across the surface, 

(t(0) = {(Jrr) = {Trr{r = a+) — Trr{r — a^))- 



(2) 



(3) 



Here 1 is the unit matrix and (• • ■ ) denotes time average 
over an optical period. Since fi = I everywhere, the 
magnetic terms of the Maxwell tensor do not contribute 
to the surface force and may be omitted. 

We introduce complex fields so that £ — Re{Ee''^*}, 
etc., where E, D, B and H are complex field vectors (de- 
noted as upright, as opposed to calligraphic typeface). 
For field components Xi and Xj, we have (XiXj) = 
^Re{XiX*}, and in particular {X^) = ^\X,\^. 

The force density on a sphere is found from Maxwell's 
stress tensor. In terms of the (complex) EM fields just 
inside the spherical surface it may be writen [17| 



ai9) 



{n'-l)(n'\Eli'\' + \E^\' + \Ej;\'). (4) 



where n = ni/n2. Superscript w indicates that the field 
is evaluated inside the droplet. We shall be working 
with circularly polarised light so that our system is ax- 
ially symmetric throughout. Projecting onto a basis of 
spherical harmonics, we may write in this case. 



aie) 



^a,P/(cos0). 



(5) 



It was found in Ref. jTTi] that the I = term is zero due to 
mass conservation and that the I = 1 term corresponds 
to uniform movement of the entire droplet and can be 
ignored, so only I > 2 are of interest to us. The orthog- 
onality relation for Legendre polynomials gives [iTj 



ai = i(2Z + 1) / de sin 9(7 (9) Pi [cost 



(6) 



We define here for future reference the quantity a 
which is the number of wavelengths (in the external 
medium) per circumference, a key parameter. 



27ra n2i-oa 
A2 c 



(7) 



(fc2: wave number in outer medium). 

In this paper we have assumed all fluid to have zero 
absorption of light for simplicity. Before going on to the 
fluid mechanics of the droplet, let us briefly discuss the 
ramifications of this assumption. Two notable physical 
effects would manifest themselves were we to include a 
small absorption coefficient for the fluids (equivalent to 
letting refractive index n have a small positive imagi- 
nary part). First, the net force on the droplet would 
acquire an addition in the direction of propagation due 
to absorption of photons, each carrying a momentum 
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hjjj. More interestingly for us, however, is the photoa- 
coustic effect resuhing from the increase in temperature 
which would result from absorbed energy, causing local 
thermal expansion of the fluid. The effect has many ap- 
plications in metrology (4^] . When the laser is switched 
on the local temperature and pressure will rise propor- 
tional to time t and absorption coefficient aabs- After 
a while, local thermal equilibrium is reached - heat is 
transported away at the same rate as it is produced - 
and the temperature rise and expansion stops. Typical 
temperature rise and fall times for micrometer sized liq- 
uid systems are about 5 — 10/iS [H, |43|. The pressure 
change creates an acoustical wave which propagates out 
of the system. After the rise time, thermal expansion, 
just like electrostriction, will contribute only to the pres- 
sure distribution, and not to the motion of the fluid. It 
can thus have a bearing on the fluid surface dynamics 
only when light is modulated on a timescale similar to 
the rise time, or short enough for acoustics to play a 
role. 

2.B. Fluid mechanics 

The optical force acting on the interface between the 
two immiscible fluids will deform the surface. The de- 
formation can be described transiently from the laser is 
switched on j^, [s^, [111 , but in the present case we shall 
assume static conditions. We assume throughout that 
the surface deformation is small enough so that we may 
keep only linear order in the deformation amplitude. On 
the other hand we will calculate the optical force as if 
the droplet/bubble were spherical, which is a reason- 
able approximation when the dielectric contrast is low, 
n2 — ni <C rii , but corrections of linear order in ampli- 
tude are expected to be of greater importance for greater 
contrast, such as for a water droplet in air as considered 
in Refs. [13, [H, HI]. (Note that these references con- 
sider the case of an optical pulse short enough so that 
the particle remains spherical throughout its duration). 
A full theory for the static deformation of trapped mi- 
crodroplets and -bubbles is under development. 

We write the surface deformation of the droplet also 
in terms of Legendre polynomials: 



(8) 



1=2 



In Ref . |T3| was considered the transient case hi {t) where 
a laser beam was switched on at t = and off again 
after a time to- The deformation coefficients in the cur- 
rent, static case is then given from the solution found in 
Eq. (10) of [l3] upon taking < — > oo while keeping < < to 
(a long time after switch-on, but before switch-off). The 
static solution is 



hi{oo) 



1 



7 P + l-2 
where 7 is the surface tension coefficient. 



(9) 



More generally, a time t > after the beam is switched 
on, the surface deformation is described for the under- 



damped and overdamped case, respectively, by [17 1 
hiit) 



hi{oo) \7/ 

hijt) 
hi{oo) 



1 - sin -fit + cos jitje w; > 



(10a) 



sinhT/t + coshr;t le"'"* 



w/ < ^ll 
(10b) 



for I > 2, where Ap — \p2 — pi \ is the difference in mass 
density and 
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' Apa^ 



[P + l- 2) 



1 I M2 (/+l)(i+2) 
pi i+1 



(11a) 
(lib) 

(11c) 



The expression for pi generalizes that given in [17|, |33j , 
which was valid only when the viscosity and density of 
the outer medium can be neglected compared to the in- 
ner. Details of its derivation are found in Appendix IB] 

The fluid dynamic theory presented is valid to lin- 
ear order in h(ff) and h!(ff). The most important cor- 
rection to this is likely to be the onset of nonlinear 
terms from surface tension, since this terms contains lin- 
earlised factors ^^1 -f [ft-'(6')/a]^, giving third order terms 
'-^ /i[/i']^/a^. The quadratic term ignored upon linearis- 
ing the Navier-Stokes equation for fluid motion, due to 
internal convection in the droplet, is of gradient form, 
and may be absorbed in the pressure, being fully ac- 
counted for [sl]. In the plots of fluid shape in the sec- 
tions to come, the maximum deformations are approx- 
imately |ft.|/a, \h!\ja ~ 0.2 — 0.5, which is pushing the 
boundaries of the linear regime a little. What is accom- 
plished is a clear presentation of the qualitative shape, 
which is the intention of the figures in question, although 
the actual shapes seen in the figure should be under- 
stood as semi-quantitative. Using amplitudes exceeding 
the linear regime for illustration purposes is standard in 
the literature, c.f. e.g. [s^ |33|. 

We shall treat both droplets and bubbles as incom- 
pressible in the following. For the latter case, this might 
require comment. For steady flow the condition that 
compressibility effects be negligible is that the fluid ve- 
locity is much smaller than the speed of sound (about 
340m/s in dry air). For the present case, fluid velocities 
are of order v ^ au: ^ Apa for underdamped case 
when velocity is the greatest. Inserting numbers for air 
surrounded by water, / = 2, velocities in the order Im/s 
are obtained, hence a Mach number of order 10~^. For 
unsteady (transient) flow such as here, a second crite- 
rion must be fulfilled (c.f. [l!] §10), that the timescale r 
of the motion must be far greater than L/cs, L being a 
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typical lengthscale and c— s the speed of sound. Here, it 
would imply aco/cz ^ 1 , again satisfied by two to three 
orders of magnitude for a micron sized droplet. 

Alternatively, we may estimate the maximum possi- 
ble change in density by regarding the maximum optical 
force density acting on the surface, which from Eq. Q 
is ~ n2(n^ — l)/o/2c ~ Sp, the change in internal pres- 
sure. Using the ideal gas law this would correspond to 
a change in density of 5p ~ 5p/R^{^T where the specific 
gas constant for air is i?air = 287J kg~^K~^, and T is 
temperature. With values used for the bubble case we 
obtain (5pmax ~ — 0.07kg/m'^ which is tiny compared to 
atmospheric air density of around 1.2kg/m^. Hence we 
may safely ignore effects of compressibility for gas-filled 
bubbles, as well as droplets. 

2.C. Overdamped or underdamped oscillations? 

Oscillations are underdamped whenever uj^ > and 
vice versa. We consider three different cases separately: 
droplet in air, air bubble in liquid, and a liquid-liquid 
system such as the oil emulsion in water. 

2.C.I. Liquid droplet in gas 

For a droplet in air, /i2 <C /i2 and P2 <C pi, so we may 
safely replace air with vacuum. The criterium for under- 
damped oscillation for any mode whose index / is then 



(2/2 



I -I) 



I > 2, 



(12) 



where 



i.e., when 



I <j^[a + 2^36 + 24A + A^ 

6V3AVM8 + UlA + 136^2 ^ 4^3 



X cos 



■ arctan ■ 



3—" 216 + 54y4 + 36/12 + A3 
A/ 4, when A > 1 
1, when A < 1 



} 



(13) 



where the last line is asymptotes for large and small A. 
When A < 25/8 = 3.125, aU modes / > 2 are over- 
damped. This is illustrated in figure [2j 

In the lingo of dimensionless numbers in fluid mechan- 
ics, 



Ai 



Gai 



(14) 



where Ga, the Galilei number, is the ratio of gravita- 
tional to viscous forces, and Eo, the Eotvos number 
(virtually identical to the Bond number) is the ratio of 
bouyancy forces to capillary forces. 



Gai — — 2" I 



Eo = 



Apa^g 



Shaded: 


Underdamped 


Unshaded: 


Overdamped 



(oil/water) 



Oil-WaterXy X 



Droplet 



. Bubble 




(water/air) 



Fig. 2. Critical damping curves for the three cases: (1) a 
droplet in air, (2) an air bubble in liqid, and (3) the two-fluid 
example of oil-in water where ^12/^1 ~ 1.2 and P2/P1 ~ 0.32 
were used. The abcissa is Ai for cases (1) and (3), and A2 
for case (2). Shaded areas denote underdamped modes. The 
interval of I values that contribute signiflcantly to the force 
density in Eq. ©, 2 < i < q, are indicated (slanted and hor- 
izontal dashed lines) for air-water and oil-water, respectively, 
when Ao = 1064nm. 



Bearing in mind the conclusion from [17|, |33| that the 
number of significant terms in the I sum is no more than 
order a (for the static case it can be significantly less), 
the behaviour of the droplet is sure to be essentially 
underdamped when a < A/ A, assuming A^ 1. 

For the example of a water droplet in air (Ap = 
997kg/m3, 7 = 0.073N/m, j/ ^ 1.01 x 10"'^m2/s), the 
value of A is 



^watcr— air 



72.5 • a/1 fim. 



For radii of a few microns, all important oscillations are 
thus underdamped. The condition that all significant 
modes be underdamped, a A/A can be written 



Ao > 



87r7i2Apz^2 

7 



this 



(Aq: laser wavelength in vacuum). For water-air, 
implies underdamped oscillations whenever 

Aq > 350 nm 

which includes the full visible spectrum. 
2.C.2. Gas bubble in liquid 

In the opposite case, ^2 ^ ^J■2 a-nd P2 ^ Pii the condi- 
tion for mode I to be underdamped is 



A2 > 



(2/ + 1)2(^ + 1)4(^ + 2) 



I > 2, 



(15) 



where 



A. = 



07 
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As shown in Fig. [2l the underdamped region makes 
up a smaller portion of the I-A2 space than was the case 
in the l-Ai space for the droplet. Clearly, viscosity has a 
somewhat greater influence on oscillations on the surface 
of bubble than on droplet. 

One may ascertain from Eq. (jlSp that all modes for a 
bubble are overdamped when 

A2 < 82.23. 

For an air bubble in water, A2 — 72.5a/l/im, hence it 
follows that there exists a range of very small radii for 
which all surface wave modes of an air bubble in water 
are overdamped, whereas those of a droplet of the same 
size are underdamped. For sufficiently large bubbles, 
however, A > 2025/8 — 253.125, all relevant modes are 
underdamped whenever Ao ^ 350 nm, just as was found 
always to be the case for droplets. Also this is indicated 
in figure [2l 

2.C.3. Two fluids 

Since droplet and bubble denote two limits of the general 
expression for the viscosity coefficient ni , corresponding 
to 1^2/ fJ-i, P2/ Pi — > and — )■ 00, respectively, it is only to 
be expected that the critical damping curve of any two- 
fluid system where the properties of both fluids must 
be included will fall somewhere in between those of the 
limiting cases. Numerically this is found to be the case. 

Due to the extremely low interfacial tension 7 for the 
oil-emulsion system, however, the nature of the criti- 
cal damping curve for this set-up is found to be of no 
relevance. Whereas underdamped modes can only ever 
occur for values of Ai in the order of 10^ or higher (see 
Fig. [5]), the numerical value found with the parameters 
used herein, chosen for ultra-low surface tension, is 



found to be 



Al 



oil-water 



0.00050 • a/1 /im 



(Ap = 320kg/m3,_y « 5 x 10-^N/m, ^oii = 1.2Mwator - 
see Refs. [l^, [2^ Hg]). The droplet radius required for 
Al ~ 10^ is thus far into the macroscopic regime. We 
can safely conclude that all surface wave modes of the 
oil emulsion are underdamped. 

In conclusion, surface oscillations of water droplets in 
air are entirely underdamped, of air bubbles can be ei- 
ther under- or overdamped, and of oil-emulsions in water 
are entirely overdamped. This holds for all microscopic 
particles illuminated by laser light in the IR, visible or 
near-UV spectrum. 

2.D. Gaussian beam profile 

The electromagnetic fields in a Gaussian beam have no 
simple closed expression, but have been found to high 
order in expansion parameter A/wp {wo is waist beam 
width) when this is small [4^. Taking the example of 
the linearly polarised fields given in Ref. jlOl (we use the 
formalism of Ref. [sij ) and adapting to circular polar- 
isation [Ey = iEx, etc.) leading order expressions are 



2Q 



El = iEl 



E. 



W 



with 



i + 2z/kwQ ' 



(16a) 
(16b) 

(16c) 
(17) 



This is a solution to Maxwell's equations, modulo higher 
order corrections. Superscript i on E and H denotes 
incident fields. This now describes the incoming fields 
from the laser beam. 

Note that Eq is not the complex amplitude of the E- 
field in the above, but is chosen so that (£^) = Eq at 
X = y = z = 0. 

We define another dimensionless parameter k, being 
the number of wavelengths per beam width 



K — kWQ. 



(18) 



In typical situations both a and k can far exceed unity. 
The field expressions of Eq. [12] are truncated at leading 
order in 1/k, hence valid for k ^ 1. This ensures that 
field expressions remain good many wavelengths away 
from the beam focus. We will see in the following that 
the condition for the approximate field expressions to 
hold within the entire droplet or bubble is 

> a, 

which is thus the criterion for the accuracy of the theory 
in the following. This is a restriction on the relative mag- 
nitudes of our three length-scales: radius, beam width 
and wavelength, 

Ao ^ „ Wo 
— < 27rn2 — . 
wq a 

2.E. Internal field and explicit force density 

In spherical coordinates the internal fields inside the par- 
ticle, at coordinate r, 6, are given by (53 | 

^'r ^^J2Y1 + l)cin^Mfiag)Yi^{n) (19a) 



l—l m— — l 
I 



l — l m— — l 



dim , ,- -.yim^S^) 

mipi{nag) — — — 

712 sm 



(19b) 



^ 1 = 1 m=-l 

- ^^■il}i{nag)d0Yim{^) 
n2 



(19c) 
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where = {9,(j)), g = rja, and the coefficients as defined 
by Barton et al. are |52| (we will define more suitable 
versions below) 

cim =^iAi„i[n'^^i{na)^^^^' (a) - f^V^(^a)C^^]^^ (20a) 
d,™ ^iBi„,[Mna)£,l^^'{a) - ni;'i{na)^l\\ (20b) 

The incident field is contained in the quantities 



Aim ^ 1 f El/Eo 

Bin, l{l + l)Ma)J Hl/Ho 



Y:^{n)dn, (21) 



where incident fields are evaluated at r = a"*" and the 
integral is over all solid angles. Here Hq = Eo/{cfio) = 

For a circularly polarized Gaussian beam the radial 
incident field is (as always truncating to leading order 
in X/wq ^ 1/k) 

El _ El sin^e*-^ + El cos 6* 
Eo Eq 



smWe 



Ex 
En 



K am Ue 



2ia 



■ cos ( 



{k? + 2iacos( 



■ exp 



ian^ cos 9 — a^{l + cos^ 



2ia cos 9 



(22) 



We see immediately that the criterion for the correction 
terms in the approximate expressions for the incoming 
fields to be small is that ^ a as stated above. 

To calculate Aim, Bim using Eq. ([2T|) we use the rela- 
tion 



Yimm 



121 + 1 {I -my. 
4tt {l + my. 



P;"(cos6l)e 



and since E^ = e^'^f{9; a, k) only m = 1 coefficients are 
nonzero: 



27r 



We hence obtain the expressions 



Aim — 



y/(2l + l)TrSml 



■$,(a,K); 



Bim = - itnA, 



hn 7 



(23a) 
(23b) 



and 



$i(a, k) = 



1 , k^Vi^PHu) 

— 77- 7^5 — 



X exp 



ian^u — q;^(1 + u^) 



2iau 



(24) 



which must be evaluated numerically. In the limit where 
the laser waist tends to infinity we find 



2i'+i 

lini $/(a, k) = + (25) 

K^oo a r—T. 1— -I 

which gives back the known plane wave case [17|, |32| . 

With this we are ready to insert our expressions to 
obtain the force density of Eq. ([4]) and project onto the 
basis of Legendre polynomials using ([6]). With some 
straightforward but tedious algebra we obtain 



e^nlEl{n^^l) 
32 



oo oo 

I c„tC* -0m (na)-0« {na)Ii 



Imn H~ ^ 



la 



Cmd*n^mina)^nina) - c^mC*-^/'m(f^a)'0^l(?^a) ^;mn|- 

with the constant coefficients 

(2; + l)(2m+l)(2n+l) 



Ml, 



Nlr^ 



m{m + l)n(n + 1) 
(2Z + l)(2TO + l)(2n+l) 



AuPi{u)Pl{u)P^{u), 



[m{m + l)n(n + 1)]^ 

(2; + l)(2m+l)(2n+l) /"^ 
[m{m + l)n{n + 1)]^ 



duPi{u) {l-u')P^{u)P^'{u) + 



Pl{u)Pl{u) 



duPl{u)Pl{u)Pl{u) 



(26) 

(27a) 
(27b) 
(27c) 



r 



and where we define more handy c and d coefficients, 
^i{a,n) 1 



Exact expressions for /;,„„ , Mimn and Nimn are given in 
appendix |X] 



ci 



Ma) n%l}i{na)i\^^' {a) - i;[{na)i[^^ 
^i{a,n) 1 



(28a) 
(28b) 
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Calculating ai is the major numerical task, whereupon 
insertion into (O or ((TU]) quickly yields the static and 
transient surface deformations. 

2.F. Note on intensity and power 

Experimental values are normally not given for Eq, but 
rather for intensity or laser power. The intensity is given 
by the Poynting vector ((•••) denotes average over one 
optical period as usual), 



/ =(5,) = {Sx H), = iRe{E x H*}, = 



n2\Ex\ 



where we have made the observation that RejiQ} = 

The power is given by the integral of the Poynting 
vector across the beam cross section, and is independent 
of z as it should: 



P = jTTWo/o = \TTeQcn2wlEl 



(29) 



Equation ([26]) can then be written 
n2h{n^ - 1) 



32c 
n2P{n^ - 1) 







rn—l n— 1 

oo oo 



(30) 



■m— 1 n— 1 



The effective power is the intensity of the laser light 
integrated over the central cross section of the sphere, 

Pc« = ^(l-e--^/-o). (31) 
In the case wq ^ a we obtain the obvious limit 

P™=^a'^o. (32) 

When comparing the case of a narrow laser beam k < a 
to one of infinite width, we will quote the effective power, 
which is the comparable unit between the two. 

2.G. Necessary numerics 

The numerical calculation involved is thus twofold. The 
coefficients Iimn , ^imn and Nimn niay be calculated once 
and for all and is numerically the heaviest part, and 
secondly the integrals $i(a, k) must be calculated for a 
given set of parameters. 

The sums in Eq. ([^51) can each be truncated after 0{a) 
terms, and similarly the number of cr/s of significant mag- 
nitude. Explicit calculated values of the coefficients of 
Eq. p7| (quoted in appendix) can be written as triple 
sums with a total of Imn terms. We see that the problem 
initially has a terrible scaling with increasing a; the cal- 
culational cost increases as a®. It is not untypical with 



values of a in the order of hundreds, and calculation 
time quickly becomes forbidding. 

A scaling closer to can be obtained by approximat- 
ing the integrals (|27p using a trigonometric approxima- 
tion for the legendre polynomials [53| . Upon intergration 
this approximation is only accurate to about 20%, how- 
ever, so the search for a quicker algorithm continues, for 
example by taking the approximation to higher order. 

3. Numerical examples 

In the following we will provide a number of numerical 
examples of the above theory. We focus on three differ- 
ent physical situations; a water droplet in air l3l| . an 
air bubble in water [36| . and an oil emulsion droplet in 
water 

3. A. Internal field intensity 

Maxwell's stress tensor is intimately connected to the 
electric field intensity, as is immediately clear from 
Eq The higher the intensity, that is, the mean 

square field, near the spherical boundary, the larger the 
local optical force density. It is of interest therefore to 
compare the electric field intensity within droplets of the 
two different dielectric contrasts found in the water-air 
system (higher contrast, n2 = n = 1.33) and oil emul- 
sion system (low contrast ni = 1.41, n = 1.06). Re- 
sults are shown in Fig. |3] for an intermediate size droplet 
(a = 15.7) chosen in respect to the wavelength so that 
scattering patterns be neither trivial nor unneccessarily 
complex. 

From the prefactor — 1 in Eq. ^ it seems that 
the local force density on the droplet surface should be 
proportional to the dielectric contrast An, but in truth 
the scaling is stronger than this because of focussing ef- 
fects from refraction. A droplet with higher dielectric 
contrast such as water in air will act as a lens, focussing 
the incoming light onto its rear interface, creating spots 
of significantly higher intensity than average. The dif- 
ference between the two cases is clear in figure [31 both 
for a beam that is narrower than the droplet 's diameter 
(Fig. [3] a and b), and for a wide beam (Fig. [3|c and d). 
In the figure we compare a situation of low dielectric 
contrast (oil emulsion in water, ni = 1.33, n2 = 1.41) in 
a and c to the higher contrast system of water in air in 
b and d. 

Because the deforming force goes like light intensity 
at the dielectric boundary, the shapes obtained can be 
qualitatively different due to the different contrast in 
otherwise comparable circumstances. A corollary of this 
is that to correctly treat systems with even higher di- 
electric contrast in the static case where the droplet is 
still illuminated after it is deformed (as opposed to a 
short pulse) , a more general scattering theory should be 
drawn upon for full accuracy. Such a theory does not 
exist to date to our knowledge, and will be the topic of 
future studies. 
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(a) Oil emulsion in water 



(b) Water in air 




(c) Oil emulsion in water 
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Fig. 3. Absolute square of electric field |E*| inside droplets with lower and higher dielectric contrast. The droplet radius in this 
example is a = Beam waist at focus in (a) and (b) is 1.5/xm (k = 11.8), and infinite in (c) and (d). Vacuum wavelength is 

Ao = 1064nm, so a = 15.7 in all cases. Contours are at the same levels in (a) and (b) (15 contours between and 2.5Eq; white 
areas: > 2.5Eq), and in (c) and (d) (15 contours between and 6Eq; white areas: E^ > 6£o)- Peak values are lASE^ (a), 
4.35£'o (b), 2.07 E'i (c), and 15.75£o (d). 



3.B. Static deformations 

The static deformations obtained for the four configura- 
tions in figure [31 are shown in the corresponding panels 
of figure m The four cases illustrate well the different 
properties observed in changing beam width and for dif- 
ferent optical contrast n. The static shapes should be 
considered semi-quantitative only, since the EM fields 
used are only zeroth order. For the oil-emulsion system, 
the optical contrast is so low that scattering is modest, 
and the approximation is not expected to be bad. For 
water-air there will clearly be corrections. A redeeming 
feature is, however, that the greatest deformation of the 
surface is found on the side of the droplet facing away 
from the laser beam, onto which the light is focussed, 
whereas the irradiated side is still close to hemispheri- 
cal with a somewhat smaller radius of curvature. Hence 



we can expect the focal point of the light to be moved 
slightly towards the centre of the droplet, and the static 
deformation therefore to be slightly overestimated by 
our theory. An improvement of this particular aspect of 
the treatment is a goal for the future. 

3.B.I. Wide vs narrow beam, high vs low optical 
contrast 

We consider now the deformations that can be ob- 
tained with a narrow beam (wq = 1.5/im as before, and 
a = 2fini) and a wide beam, applied to the low and 
high optical contrast systems, oil-emulsion and water 
droplets, respectively. 

For the oil-emulsion system we take values similar to 
those used in (2^, a total laser power of 7mW in the 
narrow beam case. Not all of the light impacts on the 
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Oil emulsion, wide beam 



'x-+y-/a 
.Xil 



(d) Water in air, wide beam 



■ iS .. 




Fig. 4. Approximate static shapes for the same configurations as in figure [S] The laser power in (a) is 7mW (Pcft = 6.8 mW), 
and the effective power in (c) is 4 times that of (a). The laser intensity in (d) is 6W/im~^ (feff = 75W) and the effective power 
in (b) is 1/3 that in (a) (25 W). Droplet radius is a = 2/im in all examples. 



droplet, so the effective power is 6.8 mW. The corre- 
sponding deformation is shown in Fig. [3^. Just like the 
squared £'-field the deformation is approximately front- 
back symmetrical, and the resulting shape is reminiscent 
of a melon. As figure|3ti' shows, the laser field propagates 
through the oil droplet almost unchanged, and the re- 
sulting deformation bears witness to the intensity profile 
being approximately Gaussian both at front and rear. 

The oil emulsion does not focus the light much due to 
low optical contrast, so when the beam width is much 
wider than droplet dimensions, the light intensity ends 
up being near uniform thoughout the droplet, as seen 
in Fig. A uniform light intensity cannot deform an 
incompressible droplet, only intensity differences can. It 
is necessary therefore to increase the effective power of 
the laser in order to obtain similar field strengths and 
hence similar deformation magnitude as for the narrow 
beam. When increasing the effective power by a factor 
4, the deformation that results is that of figure |3J:. The 
local intensity maxima seen in Fig. [3j;, on the symmetry 
axis at rear wall, and in a circular band approximately 
40° from this, show up as corresponding deformations 
in Fig. I3J;, resulting in what might be termed a beetle 



shape. 
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Fig. 5. Approximate "lemon" shape of oil-emulsion with very 
narrow laser beam. We use A — 385nm as used in [25l . |26| . 
radius a — 2.5/im (a = 54.3) and waist wq = 0.4/xm (k = 
8.7) so K^/a = 1.4, thus somewhat beyond the scope of the 
theory. 

It is prudent to note at this point that deformations of 
oil-emulsion systems are highly dependent on exact opti- 
cal contrast, laser width, laser wavelength and radius. A 
slightly different refractive index difference would scatter 



11 



the light differently in Fig. [S};, say, resulting in a differ- 
ent shape in Fig. Likewise, an even narrower beam 
than that used herein (which, remember, can only be ad- 
equately described in our formalism so long as k? 3> a), 
results in a lemon shape instead. We have approximated 
this in figure [5]- although K^/a ^ 1 in the figure, the 
shape is expected to be approximately right. In order to 
have a beam significantly narrower than a without far 
exceeding the bounds of our theory we choose compro- 
mise variables Aq = 385nm (used in [2^,[2a|), « = 2.5^m, 
Wq — 0.4/im with oil-emulsion parameters otherwise the 
same. A lemon shape results. 

Since the surface tension of a water-air interface is 
much higher than that of the oil-emulsion system, cor- 
respondingly higher laser power is required for deforma- 
tions to be comparable. We choose a light intensity that 
is similar to that of Refs. [3l| and [l^, 6W/im^^, and 
which gives a clear deformation even for such a small 
droplet as a = 2/im (much smaller than the droplets 
used by Zhang and Chang [3l| where the lowest of the 
intensities used is 4W/zm~^). While for water droplets 
this is a very small size, it is the droplet size used in 
typical emulsion systems [1^ HH] , so we use it for ease 
of comparison. The wide beam is strongly focussed onto 
the rear droplet perimeter, resulting in an intensity en- 
hancement by a factor « 16. Larger droplets (i.e., larger 
a) can result in even larger enhancements - Zhang and 
Chang report enhancement factors > 200. 

The effects of the higher optical contrast in the water- 
air case are curious to note. Contrary to the low- 
contrast case, deformations are now larger for the wide 
beam than the narrow, when effective powers are the 
same. To get the similar deformation amplitudes seen 
in Fig. |3b and d, effective laser power had to be further 
increased by a factor 3 for the narrow beam. The sec- 
ond curious observation is that the deformation is qual- 
itatively almost identical for narrow and wide beam in 
this case. The reason may be gathered from the intensity 
plots in Fig. [SJd and d. Both have intensity maxima near 
the axis on the rear surface. The intensity enhancement 
is greater for the wide than the narrow beam (maximum 
« 16Eq compared to sa AEq for this example), but for 
droplet shaping this is compensated by significant field 
intensity values along the whole droplet surface. 

We note finally a feature that is to the advantage of 
our approximate theory where electric fields for a sphere 
are used in static (trapping) deformation calculations: 
the front (illuminated) face of the water droplets are 
distorted from the spherical shapes much less than the 
rear surface, so corrections due to slightly different light 
scattering patterns are much smaller than if the opposite 
were the case (this can happen e.g. with a Bessel beam 
under special circumstances, cf. [l7|, but never with a 
Gaussian beam). 

3.B.2. Bubbles vs droplets 

The apparently minor operation of interchanging the in- 
ner and outer medium of illuminated body and exterior 
can give dramatically different results. We consider in 



this section the technologically interesting case of air 
bubbles in water, and compare it to our previous inves- 
tigations of water droplets in air with the same laser 
specifications. Internal fields and deformations for cases 
which allow comparison with the droplet and emulsion 
sections are found in Fig. [Bl to be discussed forthwith. 
One immediate consequence of interchanging water 



and air is that the prefactor of Eq. , 



1, changes 



sign so that the surface force on the bubble interface is 
compressive. A simplistic prediction would be that the 
bubble would be squashed where the droplet bulged and 
vice versa. 

The results in the two cases are more radically differ- 
ent than can be accounted for by this change of sign, 
however. The reason is found in the different light scat- 
tering properties of droplets and bubbles. While the 
droplet acted as a positive (converging) lens, focussing 
the light beam onto or near to its rear boundary, the 
bubble acts as a negative (diverging) lens, spreading the 
light out and away from its interior. This is strikingly 
visible upon comparing Fig. [SJi with Fig. [3Ji. 

In figure [5] we study the case of an air bubble in water 
for comparison with the water droplet case in figures [3] 
andUl Perhaps the most striking observation is that the 
bubble is deformed much less than the droplet under the 
same conditions, which is particularly noticeable for the 
wide beam. 

In order to visualise the kind of shapes which result 
with a narrow and broad beam, we also show the cases 
where the intensity is boosted sufficiently for deforma- 
tions to be properly visible. Consistent with the obser- 
vation that light is scattered away from the rear interface 
rather than focused onto it, the bubble is more clearly 
deformed on its front surface than its rear, in contradis- 
tinction to the case for water droplets. The shapes are 
qualitatively different from the ones previously observed, 
reminiscent of a sea urchin and an acorn, respectively. 

A final note about optoacoustics is warranted. When 
considering a laser trapping a bubble in a surround- 
ing liquid, the incoming and scattered laser beam will 
cause the liquid to contract slightly in regions of higher 
laser intensity. As the laser is switched on, a near- 
cylindrical sound wave is generated originating from the 
beam, propagating outwards. A detailed discussion may 
be found in Ref. 5J|. We have previously ignored the 
effect of electrostrictive compression, because it will be 
countered by an increase in mechanical pressure on an 
acoustical timescale, and plays no role in the dynamics 
after a time which we considered to be negligibly short. 
With a surrounding fluid, a soundwave originating from 
the laser beam can be reflected off nearby surfaces and 
arrive back at the bubble system some time after its 
emission. We have still neglected this effect since the 
sound wave's energy is distributed over a cylindrical sur- 
face, so a wave that has travelled far (and thus arrives 
late) will have comparatively little energy. In an un- 
bounded fluid regime, it is of no consequence for later 
time, since the sound wave propagates away from the 
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Fig. 6. Bubble of air in water, for the same cases as considered in figures [3] and ID a narrow beam (a-c) and a wide beam (d-f). 
Field intensity contour plots a and d are scaled the same as corresponding plots in figure [3] for direct comparison. Figures b 
and e show deformations using the same laser powers as used for water droplets in figure |4]b and d, respectively. Panels c and 
f show the same deformations as panels b and e, respectively, but with power boosted by the indicated factor. 



system studied. 

3.B.3. Comparison with experiment by M0ller and 
Oddershede 

We are in a position now to compare the deformations 
calculated here with those measured by M0ller and Odd- 
ershede (35} . In their experiment the change of maximal 
radius of the droplet in the cross-beam plane, rj^, with 
increasing laser power. By using a simple theory based 
on apparently quite severe simplifications, they are able 
to predict the relation 



(^front + -Pback) 

127r7 



-AP 



where i^front.back IS the force on the front (rear) half of 
the droplet. In the system used in [35[ the constant A 
is calculated to be 3.33^mW~^, and is found to be in 
surprisingly good agreement with the theory, especially 
since a simplified expression for the optical force is used, 
with the assumption that <C a, which is reasonably 
well satisfied in their system (they have wnjy 1/Ltm and 
a = 2.3/^m to 2.9/^m). The system in Ref. |35{ thus most 



resembles that shown in Fig. [5l producing a "lemon" 
shape. 

Let us briefly compare the linear constant with the one 
calculated and measured by M0ller and Oddershede. We 
have performed the calculation of droplet shapes using 
the data from [35] (7 = 1.9/iN/m, ni — 1.366, 71.2 = 
1.330). We find that the shapes obtained, and the con- 
stant A in particular, depends quite sensitively on the 
width of the laser waist, a parameter not given precicely 
in Ref. 35|, although wq < Ifim is indicated. For laser 
waists greater than droplet radius (P interpreted as Peff 
in that case), we find A ^ 2.3/im/W, whereas numbers 
as high as A ~ 4.5/im/W are found when the beam waist 
is made so small that k} ja ^ 1. As a general note, some- 
what higher A than the experimental one are found for 
narrow beams producing lemon shapes, whereas better 
agreement is found for beams producing ellipsiod shapes. 

With its narrow- waist laser, the geometries of Ref. [11] 
lie on or outside the boundary of validity for our the- 
ory, which is that k} ^ a, hence perfect correspondence 
should not be expected. Their observed value lies in 
the range of values predicted by our theory. In effect. 
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the approximate field expressions of Eqs. ([16)) are valid 
when the beam does not broaden perceptibly compared 
to its waist within the confines of the droplet. It is likely 
that this does not hold for the experiment of M0ller and 
Oddershede. 

The present theory, however, would seem to indicate 
that the good correspondence between theory and exper- 
iment found in 35 ] might be partly down to luck. The 
shape assumed by those authors is an ellipsoid, whereas 
the shape produced by a very narrow beam when the op- 
tical contrast is low we find rather to resemble the lemon 
shape of figure [5l Moreover, three different droplet radii 
are investigated, and our results would seem to indicate 
that slightly different values of the slope A are to be 
expected for different values of K^/a. Ref. •SS^ does not 
give a quantitative bounds on the measured slope A at a 
given confidence, so quantitative comparison with their 
experiment is difficult. 

3.C. Transient deformation: trapping beam vs short 
pulse 

We finally analyse the dynamics of droplet deformations 
in laser beams. We wish to study and compare the dy- 
namics of two different set-ups analysed in previous sec- 
tions, a water droplet in air and the oil-emulsion in wa- 
ter. Dynamics of a gas bubble in water is found to not 
be very clearly and qualitatively distinguishable from 
that of a water droplet, and will not be given specific 
attention in this context. 

For a laser pulse abruptly turned on at i = and off 
again at t — to, Eqs. pUj) generalize to [iTj 



hijt) 
hi{(X)) 



=1 - f— sin7,t-f cos7;i')e"'''* 
-hcos7z(t- to)] e"^' }, uji>iii (33a) 



and 



hijt) 
hi{oo) 



(p- sinhr,t -I- coshr,iU~'''* 



1 - 



+ e(i-<o){] 

+ coshri(t-to) c"^'(*"*«'|, LUiKfii (33b) 



^sinhr,(t-<o) 
1 1 



for t > 0. 

An illustration of the dynamics of shapes is shown 
in Fig. [7^, where a series of shapes of a water droplet 
is shown. The parameters are as in Fig. |4li, but with 
somewhat higher intensity, Iq = lOW^m"^, in order 
that the shapes be more clearly visible. In the top row 
is seen the case of a static beam switched on at < = 0, 
and in the below row it is switched off again at to = 
250ns, whereupon the droplet returns to spherical shape 
when oscillations have died out, on a time scale of a few 
microseconds. 



We found in Section 12. CI that the dynamics of water 
droplet and oil-emulsion are very different. The motion 
of the former is in practice always underdamped, the 
latter overdamped. Moreover, the time scale of a given 
mode I is very different for the two. For an underdamped 
mode, a relevant timescale is its lifetime 



droplet 



1 



4.0/iS 



droplet 



{2P -l-l) 



(34) 



for water. The time scale of the overdamped oil emul- 
sion, however, follows from Eq. (1101) wherein 3> LL>f, 
and 



^emulsion 



7 lf{l) 



0.012s 



(35) 



where v is the mean kinematic viscosity and f{l) is a 
function of order unity (exact details can be derived from 
Eq. (Uni) if of interest). 

A striking difference between the water droplet and 
the oil emulsion is the timescale involved for the fluid 
dynamics. The motion of the water droplet reaches equi- 
librium within a couple of microseconds, whereas the oil 
emulsion takes up to a tenth of a second to find its fi- 
nal form. This difference is illustrated in figure [7]d and c 
where the motion of the front and rear tips of the droplet 
are plotted. Above the abcissae the rear point position 
r{d = 0)/a (where the light exits) and below the front 
r{d = 7r)/a (where the light enters). 

As noted also in Ref. [13, [12, , the higher I oscilla- 
tion modes are short-lived compared to the lowest order 
modes, as also Eq. shows clearly. In Fig. [TJi and 
e we show the time evolution of the four lowest oscilla- 
tion modes, as given by the right hand side of Eq. p3ap . 
By the time the I = 2 mode reaches its maximum, only 
/ = 3 still contributes significantly to oscillations. In- 
deed, for the majority of the oscillation period of about 
3fj,s in the static beam case, only the / = 2 mode con- 
tributes significantly to the time variation. That is not 
to say, of course, that the final shape of the droplet can 
be adequately described by the I — 2 shape perturbation 
alone. 

Note that for a pulse of 250 nanoseconds we may ig- 
nore compressibility effects due to elcctrostriction (see 
Ref. fioj). This becomes of importance when pulse du- 
ration is in the order of traversing time of a sound wave 
across the diameter, ts = 2a/ Cs where Cg is speed of 
sound. For liquid and gas bubble respectively, this time 
scale is approximately 4 and 12ns, which is much shorter 
than our pulse. 

4. Concluding remarks 

We have laid out a comprehensive theory for linear de- 
formations of fluid micro particles in a laser beam of 
Gaussian profile. Three types of fluid systems were con- 
sidered and compared: a water microdroplet in air, an 
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Fig. 7. Dynamics of deformation for a water droplet, comparing the cases of static beam and a short pulse, (a) Water droplet in 
wide beam, as in Fig. |4|i, but with slightly higher light intensity, Iq — 10W/im~^. Static beam (upper row) compared to beam 
of duration 0.25^s (lower row), (b) Time development of droplet elongation for the above shapes, illustrated by r{0 = 0,t)/a 
(rear of droplet, plotted above abcissa) and r{9 — n,t)/a (front of droplet where light enters, plotted below abcissa). Static 
beam turned on at ty = (solid) and pulse of duration to = 0.25^s (dashed) are shown. Panel (c) shows the same as (b) but 
for the overdamped case of an oil emulsion in a narrow beam, same set-up as Fig. |4^, but with power boosted to 12mW, for 
static beam switched on at t = (solid) and a pulse switched off again at to = 10ms (dashed) . Panel (d) and (e) show the time 
development of the four lowest deformation modes, as given in Eq. (|33|l . for a static beam switched on at i = (d) and a pulse 
switched off again at to = 0.25ps. The legend in panel (e) is for both panels (d) and (e). 



air microbubble in water, and an oil-emulsion in wa- 
ter. The latter system has found several experimental 
applications due to its extremely low interface tension 
coefRcient. 

The fluid dynamics of the three types of fluid sys- 
tems are surprisingly different. We show that dynamics 



of a water microdroplet in air is always underdamped, 
whereas the oil-emulsion always shows overdamped dy- 
namics since the surface tension, which acts as a "spring 
constant" in the harmonic oscillator-type equations of 
motion, is very weak compared to the friction term from 
viscosity. A microscopic bubble in water can have signif- 
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icant contributions from both underdamped and over- 
damped modes, depending on radius and laser wave- 
length. 

Although the water droplet differs from the bubble 
in water only by an interchange of materials in the in- 
ner and outer region, their deformations from radiation 
pressure upon laser illumination are strikingly different. 
All else being equal, a laser which significantly disfigures 
a water droplet will typically hardly change a bubble at 
all. The reason being that a droplet acts as a posi- 
tive lens, focussing the incoming light onto its shadow 
face, whereas a droplet is a negative lens, scattering the 
light away from its interior. A much higher laser power 
is therefore required to visibly deform a bubble than a 
droplet of the same size. 

We furthermore analyse the kind of shapes obtained 
when particles of the three different systems are illu- 
minated by a wide laser beam (beam width 3> radius) 
and fairly narrow beam (beam width 3/4 of radius), re- 
spectively (the present theory is unable to handle beams 
which are much narrower than the droplet). Qualita- 
tively different shapes are found, as we show in Figs. U 
and [HI as may be summed up as in Table [T] 





Droplet 


Oil-emulsion 


Bubble 


Wide beam 


Egg 


Beetle 


Acorn 


Narrow beam 


Egg 


Melon/lemon 


Sea urchin 



Table 1. Qualitative shapes of fluid particles under illumi- 
nation by narrow and wide laser beams. The shapes refer to 
those seen in Figs. |4l [5] and [6] 

The light scattering pattern within the droplets and 
bubbles has everywhere been calculated as though the 
particle were spherical. This is near exact for laser pulses 
which are so short that the fluid mechanical motion hap- 
pens after the light is off, but it is obvious that correc- 
tions must result in the case of a static beam, and the 
shapes calculated in this case must be considered semi- 
quantitative only. We expect this approximation to be 
fairly good in the case of oil-emulsion systems, where the 
dielectric contrast is small and focussing only very slight, 
as seen in Fig.[3]a and c. For water-air systems, however, 
the light intensity pattern is rather sensitive to exact 
shape. A redeeming feature is that the deformation of 
a water droplet is most prominent on the rear (shadow) 
side, and less on the side of the droplet impacted upon 
by the light. From the egg-like shapes in Fig. [3] b and 
d (which are exaggerated for illustration purposes) the 
correction would be due to the somewhat smaller radius 
of curvature on the side of the droplet facing the light, 
causing the focal point to move away from the rear wall 
somewhat. The expected effect might be a somewhat 
diminished deformation, yet qualitatively similar. We 
intend to return to this question more carefuilly in the 
near future. 

The dynamics of a water droplet and an oil-emulsion 
were finally analysed. While not specific to a Gaussian 



beam profile, such a comparison of dynamics does not 
exist in the literature to our knowledge, and is of obvious 
interest to laser manipulation of different kind of fluid 
systems such as those considered. Because of the very 
low interfacial tension of the oil-emulsion droplet, its dy- 
namics are slow compared to the water droplet. For the 
example of a droplet 4/im across, the relevant timescales 
are tens of milliseconds and a few microseconds, respec- 
tively, i.e., they differ by 4 orders of magnitude. We 
analyse more closely the dynamics of the water droplet, 
being underdamped. The longest wavelength oscillation 
mode I — 2 was found to be the most long lived. In- 
deed, for approximately the latter 2/3 of its lifetime, 
this mode is the only one contributing significantly to 
dynamics. The static shapes of the droplets, however, 
can have significant contributions from modes all the 
way up to / = a (a is the number of wavelengths per 
circumference), although surface tension tends to smear 
out the short- wavelength perturbations and lending em- 
phasis to lower I modes. 
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Appendix A: Exact values of I,M,N coefficients 

The coefficients in equation (|27|) can be evaluated ex- 
phcitly. Let A^imn and Afimn be the integrals in 

equation pT]). We find the expressions 

(X,A^,AA)z„„= § (z,m,n)ir^"' (Al) 

/'m'n' 

where we use shorthand 

J„,=EE E0(™)(:,) 

/ l+l' -1 \ / m+m' — 1 \ / n+n' -1 \ 

H • )( .» )( n ) 

X e{l + l')e{m + m')e{n + n'), (A2) 
A^l + m + n, A' = I' + rn' + n', (A3) 

and use the even/odd selectors 

I 1, n is even , , | 0, n is even , . ,, 
10, n IS odd II, n IS odd 
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The specific summands are 



•I'm'n' r)A+2 



e(A)(A'^ 



TO, 



I'm'n' rjA+l 



Imn 



I) 



m n 



(A5a) 

2m'n' — m' — n' — 1 



A' + l 

(m'-l)(n'-l)i 



A' - 1 



A' - 3 



(A5b) 
(A5c) 



The even/odd prefactors removes half the terms. An- 
other observation is that all terms are zero for which 
one out of I' , m', n' is greater than the sum of the other 
two. 

Explicit evaluation of these sums is numerically rea- 
sonable for moderate values of to, n, yet since the num- 
ber of terms is of order Imn, and and since the sum in 
Eq. (I26p requires terms up to / ^ a and m,n I, the 
full evaluation of cr; requires a number of terms which 
scales as a^. This quickly becomes forbidding for large 
radius-to- wavelength ratios. This is the geometrical op- 
tics limit, however, where other, approximate methods 
may be employed. 

Appendix B: Viscous coefficient 

Following the recipe of Brevik and Kluge [s^l we derive 
the approximate expression for the viscous damping co- 
efficient, Eq. Illcl Since the velocities are small, we can 
make use of a velocity potential as is standard in linear 
gravity wave theory, v = V(/3. Incompressibility dictates 
V^(^ = 0, and hence the internal and external solutions 
can be expanded in modes. Then 



(Pint 



:^GPr(cos0)e™*^)' 



1=0 



1=0 



The interface is characterized by the equation S = a -\- 
h{6) — r = 0. The kinematic boundary condition at the 
interface is then given by [ssj 

— = 5+(v.V)5 = 

which, when keeping only terms of leading order in the 
small amplitude h, gives simply 

'f\r=a=Vr\r=a (Bl) 

wherewith we match the inner and outer solutions to 
find 

I 



D, 



l + l 



Following Ref. [33| we work in analogy to a damped 
harmonic oscillator. The time average of complete me- 
chanical energy is given as i?mech = / d'^xpw^ ([48] §25), 
that is, 



^^moch = Pi / d a;(V((5i„t)2 + p2 / d x(V(/3ext)^ 
J int 



Consider in detail the internal part: 

E 



Enlch =^Pl«EE 



I 



1=0 l'=0 m=-l m' = -V 
1 r2TT 



X / Agg 

'0 



l+l' 



Corrections to the integrals due to surface deformations 
(upper limit of g integral is generally 1 + h/a) give rise 
to second order terms in the differential equation and 
are neglected. Using the integral relations quoted and 
derived in appendix [C] all non-diagonal terms are found 
to be zero, and the result is a simple sum over / and m 
only: 



E 



mcch 



oo / 



E 21 + 1 ii-my: 

1=0 m=-l ^ ' 



A completely analogous calculation gives the external 
mechanical energy, and the sum is 



E 



OO / 



mcch 



. . (? + to)! / 

1=0 m=-l ^ ' 



IP2 
l+l 



There is thus a mechanical energy uniquely associated 
with coefficient Cim, hence with mode l,m. 

The time average of the dissipated energy within a 
volume may be expressed in terms of surface integrals 

bounding said volume, i?mech — — § d^x^Vv"^ ■ n ( [isj 
§16) where n is a unit vector normal to and pointing out 
of the bounding surface. There is no contribution from 
the surface at infinity, so we obtain 



E 



mcch 



I 



Stt 



\Cl.n?Yl E ' 



{I + to)! 



1=0 ' 



{I -my. 



+ /i2(/ + l)(/ + 2)]. 



Here / dil is an integral over all solid angles. Noting the 
relation ([HI §25) 



fJ-ln 



E- 



mcch,/m 



-£^mcch,/r; 



(B3) 



equation (|llc[) results. As for the special case consid- 
ered by Brevik and Kluge [ssi], the damping coefficient 
does not depend on m (only modes to = 1 are excited 
by the Gaussian, circularly polarized laser anyway, but 
the general result will be useful for future, more general 
work. For the linearly polarised laser considered in [s^ , 
modes to = ±1 are excited, whereas an order toq Bessel 
beam excites modes to = TOq [T7|). 
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Appendix C: Some integral relations 

To derive Eq. (|B2p we require the well known relations 



d0e^(™-™')^^27r5w 



5w. 



In addition we require the integral 



Axf{x) 



rPr(x)Pr(x) 



Using Legendre's equation 



(CI) 



(l - x^)dl - 2xda; + l{l + I) - 



Prix) = 



and 



{d^prmpn = duprpn-pndipn-pi^id'.pn. 

we may write 

/(^) = l{dAi~x')d,prpr+[i{i+i)+i'ii'+i}]prpn- 

The integral of the first term from —1 to 1 is clearly 
zero from the fundamental theorem, and so, inserting 
into Eq. ([CT]) we find 



2/(^ + 1) (/ + m)! 
2/ + 1 {I -my. 



Siu. 



(C2) 



This integral was worked out by R. Kluge and used, but 
not quoted, in Ref. [33| (I. Brevik, private communica- 
tion). 
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